Mesoscopic Luttinger Liquid Theory in an Aharonov-Bohm Ring 



(N 
O 

o 

(N 

o 
O 



o 

(N 
> 

O 



O 



o 



X 



Mun Dae Kim^'-^, Sam Young Cho^^'^, Chul Koo Kim^'^, and Kyun Nahni-^ 
^Institute of Physics and Applied Physics, Yonsei University, Seoul 120-749, Korea 
^Center for Strongly Correlated Materials Research, 
Seoul National University, Seoul 151-742, Korea 
^Department of Physics, Yonsei University, Wonju 220-710, Korea 
(Dated: February 1, 2008) 

A careful study on the mesoscopic persistent current in a Luttinger liquid ring is carried out. It 
is shown that discreteness plays an important role in calculating the persistent current caused by 
the magnetic flux. At zero temperature, the current is shown to be independent of the interaction 
even when g = (172 — 94)/2 is not zero. The current becomes enhanced at finite temperatures with 
respect to the non-interacting case, when the parameter g is positive. 



With rapid development of nano-fabrication technol- 
ogy much attention has been paid to the mesoscopic 
system which exhibits strong quantum phase coherence. 
Many efforts have been devoted to incorporate the inter- 
action effect between electrons because they are strongly 
correlated in low dimensions. The Luttinger liquid(LL) 
theory ||, ||, ||, ^ is a standard method to investigate 
ID interacting electron systems. Therefore, it became a 
starting point to study the mesoscopic property such as 
the persistent current(PC)[||. 

Here, we consider the zero mode current in a meso- 
scopic LL ring without the Aharonov-Bohm (AB) flux and 
extend it to the persistent current produced by the AB 
flux. In the continuum limit of LL theory, the current is 
calculated through the continuity equation. In this ap- 
proach, the mean current carries a prefactor, vj, which 
is renormalized by interactions!^. However, for a meso- 
scopic ring, the momentum level is discrete and, thus, 
the continuum approach is not valid. The mean current 
must be obtained from the zero mode. We show that 
the zero mode current operator does not depend on the 
interaction parameter. 

In this study, we consider the case when left- 
right symmetry is absent and show the continuum 
approximation Q is not appropriate for the PC in the 
Luttinger liquid model, where discreteness plays a cru- 
cial role. In the present calculation, we take discreteness 
into the calculation explicitly and show that the calcu- 
lation yields new results that, at zero temperature, the 
expectation value of the PC does not depend on the LL 
interaction with g > 0. At finite temperatures the in- 
teraction enhances the amplitude of PC with respect to 
that of the non-interaction case. 

We, first, study the current in a mesoscopic LL ring 
without AB flux. Here and after, we consider the spin- 
less particles, because the spin degrees do not give any 
qualitatively different results for the amplitude of the cur- 
rent. The current can be defined through the continuity 
equation such that ^p{x) = i[H,p{x)] — Vj{x), where 
p{x) = J2p^^^'^ pip) = ^j^dxifix) is the bosonic density 



operator in the LL. Then j{x) is represented by 

j{x)^t[H,J2-e'''^Pp]+£. (1) 
■^-^ in 
p 

where c is an operator independent of x. The mean cur- 
rent operator, /, can be defined by / = dx' j{x'). 
In the limit when the level spacing Ap — > 0, the mean 
current / is given by / = limp_,o '^[H, Pp], which results 
that / = vjj^ with the interaction dependent prefactor 
vj. However, for a mesoscopic wire, we cannot assume 



the continuum limit and, 
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thus, should sum over the dis- 
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since [H, N] — 0. Here, / = c is the uniform zero mode 
current which can not be obtained through the continuity 
equation, but only through the definition in terms of total 
momentum of electrons. 

When the total particle number No is odd, there is 
no current excitation at the ground state. The num- 
ber of particles at the right(left) branch, No+{No-), 
becomes No+ = No- = {No — l)/2, where we sub- 
tracted 1, the number of particle at the zero momen- 
tum state. For even No, No+ — {No — 1 + kj)/2 and 
No- = {No ~ I ~ kj)/2. If the right(left) branch has 
one more particle at the ground state, kj = 1(— 1). The 
case when No is odd, corresponds to kj — 0. At fi- 
nite temperatures, current excitation number iV-|_(iV_) is 
added to No+{No-). For the non-interacting case, the 
zero mode current, /, can be calculated in terms of total 
momentum, —j;J2n^-^N'l-N^ T"""' represented by 
the excitation numbers, 
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(3) 



Here, J = N^ — iV_ is the current excitation number and 
N = A^_|_ + 7V_ the charge excitation number. For an in- 
teracting case, this range of summation remains correct. 
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if internal interactions do not change the total momen- 
tum of the system. The zero mode of the Hamiltonian 
exactly satisfies this condition. The interaction effect is 
included in obtaining the expectation value, (/). How- 
ever, the current operator should have the interaction in- 
dependent fermi velocity, vp, as a prefactor instead of vj. 
Here we define the fermi velocity such that vp = "^x"' 
since this corresponds to the product of the particle num- 
ber at one branch with the momentum discreteness. 

In the LL theory, the non-zero and the zero mode in the 
Hamiltonian are decoupled and the zero mode is given by 
^ {VP -f ^)( J-f + ^ {VP +^-^)N\ The contri- 
bution of the flux in the kinetic part comes from changes 
in the momenta of electrons in accordance to the twisted 
boundary condition and changes continuously. But, in 
the interaction part, excitation energy is determined by 
the current and charge excitation number which are inte- 
gers. In the continuum field theory, the particle distribu- 
tion is continuous and, thus, the interaction energy can 
also be continuous. However, since we are now studying 
the mesoscopic regime, where the particle discreteness is 
crucial. We show in the following that a small level shift 
without charge or current excitation does not cause any 
change in the interaction energy. 

In order to account the level discreteness correctly, the 
twisted boundary condition should be implemented from 
the beginning of the bosonization process. Also, a care- 
ful analysis on the discreteness of the particle distribu- 
tion is required, since the parity effect comes from the 
discreteness of the particle numbers. When the fcrmion 

field is expanded, '0r(a;) = (x)^ Sfel-oo '^"^'''^Cfcr, the 
twisted boundary condition gives a condition that k = 
^{n + (^/(/)o)0. Here, n is an integer. With this 
modification on k, one can proceed the bosonization pro- 
cess exactly same as when = 0. Considering the flux 
range, — ^ < j-^ < X'^^ note that positive magnetic 
flux causes an upward(downward) shift of the right (left) 
branch electrons. Now, the zero mode contribution to 
the kinetic energy excitation can be obtained from the 

prescription = J2r ^kr "^^^^ ^kr = ^ ■ cl^Ckr 
where :: signifles the normal ordering. 

To obtain the zero mode bosonic form in H^, we con- 
sider the ladder operator which increases the number of 
electrons above the Fermi level, such that Ur\Nr, N-r) — 
\Nr + l, N-r)- The ladder operator can be constructed as 
Efc clpS [rk - {kp + fiNr- ^))] 111 I, where Nr repre- 
sents the Nrth excitation state in the r branch. Here, the 
energies are measured from the reference level, vpj^. The 
twisted boundary condition shifts the momentum levels 
by and the ladder operator creates particles at the 
shifted levels such that. 



each branch is given by 



rk-[kp + ^{Nr-l+rp, 



(4) 
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For the even parity, application of the magnetic flux re- 
moves the degeneracy of the uppermost particle occupa- 
tion. Since positive (j) shifts the level to the right side, 
the uppermost particle takes a level at the left branch. 
Thus, the energy cost for the excitation in the left branch 
become increased by the level width, vp^, for each ex- 
cited particle. Considering this effect and Eq. (||), we 
obtain for the kinetic part 
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2{J + kj)'^-f\ , (6) 



where kj = — 1(-|-1) for positive (negative) for the even 
parity and for the odd parity. 

The zero mode in the interacting part, H°^^, is di- 
rectly deducible from the Hamiltonian. The wave vec- 
tor shift due to the twisted boundary condition makes 
change in the bosonic density operator as p(p) = : 
A , c, , 2x :. We, thus, can write H°„f for the 
p = mode from the Hamiltonian, 
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The zero mode part of the kinetic energy excitation at 



where the normal ordering subtracts the inflnite ground 
state density of type-r fermions. We observe that Nr = 
Erfc : 2xC^fc+JL2iL : is an integer number of 

the excited fermions for any value of 0. Also, Ny is a 
good quantum number because [i/, A^^] = and does 
not change continuously as implied in the continuum ap- 
proximation. Here, we note that the interaction param- 
eter 52(54) is introduced to describe the forward scatter- 
ing between different (same) branch particles excited from 
the fermi level. In order to perform the thermal average 
of the PC, we consider the energy levels of the excited 
states. Because the excitation energy in the interaction 
part is completely determined by the charge excitation 
number. A, and the current excitation number, J, the 
parity effect does not appear. Therefore, we obtain for 
the zero mode of the interaction part, 

Here, flux </) does not affect the interaction Hamiltonian. 
The PC is known to be determined by the current ex- 
citation number J and the AB flux 0. Since the charge 
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excitation degree is decoupled from the current part, it 
is sufficient to consider the total current excitation con- 
tribution only. 
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{j+kjr+2{j+kj)^) + ^j' 



(9) 
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where g — ^ ■ 

In the LL Hamiltonian, the zero modes are decoupled 
from the non-zero mode. Here, we write down the zero 
mode of the partition function 

(10) 

where terms irrelevant in the current calculation are ab- 
sorbed in constant C. When there is no charge excitation, 
the eigenvalue of operator, J, becomes an even integer, 
2m. It is necessary to carry out the exact summation 
process instead of the integration process of the contin- 
uum approximation to incorporate the AB flux effect and 
the parity effect correctly. The result is given by 



C^3 — p^(fcj + -f),e ^ 
\ L (ho 



(11) 



where C is a constant which does not depend on the flux 
(j) and 9^{v,q) = X]^-oo "^'^"^ Jacobi theta 

function. 

The current can be obtained as before calculating the 
total momentum, —j^ '}2,n=-^N„^-N^ T" (" ^) ' ^"^^ 
represented by the excitation numbers. 
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Here we note that lo 



eVF 
L 



(12) 



in unit h = I. With Eq. (|T^) and Eq. (|T2|), we readily 
obtain the expectation value of PC given by 
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which has a flnal form 

(I) = -h 
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where T* 
we used 



is the characteristic temperature. Here, 



^"^^(*"''^)=47rzV(-l)' 
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which is valid when \a\ < — lngpT[| 



■sinh27rna, (15) 



When the interaction preserves the left-right symme- 
try, = 54 i-e. g = 0[D, we observe that there is no con- 
tribution from the interaction at all temperatures. This 
relation, = 0, is not satisfied generally. Chiral edge 
currents in a quantum Hall annulus gives an example [p2[. 
Symmetry breaking due to an AB flux penetrating annu- 
lus induces the PCs. The inner and the outer edge cur- 
rents consist the left and the right branch respectively in 
a conventional LL. The interaction strength between the 
edges can be represented by the 172 parameter. Because 
the interedge interaction is weaker than the intraedge 
one, we observe that 5 > 0. 

In order to determine the PC when g ^ 0, it is 
necessary to determine the possible values of fcj. Us- 
ing Eq. (O) and the valid range of Eq. (15) with 
a = ^^^p(fcj + 1^) and q = exp(-2^), we find 
kj^O or -1 for < 0/0o < 1/2 and kj=0 or 1 for 
— 1/2 < (/)/0o < 0. The constraint on the current num- 
ber further reduces to kj = for the odd parity and 
kj = — 1(-|-1) for positive (negative) 4> for the even parity. 
We see that this is consistent with the values of kj which 
we obtained from the physical argument previously. 

We examine the zero temperature limit of the PC for 
positive g using Eq. (^^ and obtain results which are 
same to those in non- interacting caseS: 
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odd parity, 

even parity. (16) 



The above results clearly show that the PC is not renor- 
malized by interaction at zero temperature contrary to 
the results of previous calculations . This result is con- 
sistent with microscopic derivations |13|. 

For finite temperatures, the PC shows a qualitatively 
different behavior from the non-interacting case as shown 
in Fig. 1. We, first, consider the even parity case. It is 
sufficient to consider the positive </> regime(A;j = — 1). We 
can see from the partition function in Eq. (|l^) that, in 
the non-interacting case at the ground state, m = 0. The 
first excited state, m — +1, has a lower energy than that 
of the state, m = — 1. This reduces the value of current 
because value of the current has the opposite sign of the 
momentum as shown in Eq. (|T^. The interacting term, 
however, does not favor any particular current direction 
because of the equal energy cost for m = ±1. There- 
fore, inclusion of the interaction term makes the current 
reduction smaller and, thus, the amplitude is enhanced 
compared to the non-interacting case(Fig. 2). For the 
odd parity, the PC is just the translation of the even par- 
ity case by 0o/2 as can be seen in Eq. ( p^ ) with kj=0. 
The reduction of amplitude is also smaller than that of 
the non-interacting case. We note that the translational 
symmetry between the odd and the even parity case and 
the magnetization direction due to the PC are consistent 
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scopic ring. In mesoscopic systems like AB ring, the level 
discreteness plays an important role. We obtained the 
zero mode in the LL Hamiltonian by exact summation in- 
stead of continuum integration. At zero temperature, the 
PC generally does not depend on the interaction, which 
is consistent with microscopic derivations. At finite tem- 
peratures, amplitude of the PC is enhanced with respect 
to the non-interacting case and its behavior confirms the 
Legget's theorem. 

This work was partly supported by the Korea Re- 
search Foundation(99-005-D00011) and the Korea Sci- 
ence and Engineering Foundation through Center for 
Strongly Correlated Materials Research(SNU). 



FIG. 1: PCs at several temperatures for the even parity, when 
g = 0.05. The figure in the inset shows the odd parity case. 
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FIG. 2: PCs for several interaction parameter g values at 
T/T* = 0.5. The amplitudes of PCs are enhanced as g in- 
creases. The zero temperature current for an arbitrary g is a 
guide for eye. 



with the Legget's theorem p4|. 

Several groups have carried out numerical calculations 
for PCs in interacting systems with disorder 0. 
In those cases interactions preserve the left-right sym- 
metry and, thus, they correspond to the case, = in 
LL scheme. Lanczos methodjl^ showed that, for a clean 
ring, the amplitude of PC is independent of the inter- 
action strength away from half-filling. In the presence 
of disorder, however, the amplitude become dependent 
on the interaction strength. Lanczos method for a dirty 
ring^^ and Monte Carlo simulation for a electron inter- 
acti ng with a diffuse environment demonstrated that 



the amplitude is supressed, while it is enhanced in DMRG 
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